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Abstract. We study the interaction between dark matter and dark energy, with dark energy 
described by a scalar field having a double exponential effective potential. We discover 
conditions under which such a scalar field driven solution is a late time attractor. We observe 
a realistic cosmological evolution which consists of sequential stages of dominance of radiation, 
matter and dark energy, respectively. 
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1 Introduction 

The “late time ” accelerated expansion of the universe was discovered in 1998 [1]. A repulsive 
gravity-inducing (in other words, with negative pressure) dark energy component can cause 
such accelerated expansion. A cosmological constant, having equation of state w = £ = — 1, 
is the simplest candidate for this dark energy [2, 3]. 

Although such a model is fairly concordant with observations, it suffers from two 
major issues: smallness of the cosmological constant, and the coincidence problem. Its ef¬ 
fects are becoming just noticeable in the recent history of the universe [4]. A way to alleviate 
these problems is to explore frameworks that cause an effective cosmological constant to 
dynamically arise out of interaction between dark energy and dark matter - thereby repro¬ 
ducing the late time accelerating behavior quite akin to that produced by the cosmological 
constant [4-6]. 

The present article describes dynamical features of tracker solutions that can be ex¬ 
pected to arise from the dynamics of moduli fields that naturally occur in string theory [7]. 
These fields provide a parametrization of a compactified manifold, while their vacuum ex¬ 
pectation values determine the effective four dimensional gauge and gravitational coupling 
constants from the unification scale. It has been shown that in a universe dominated by 
other matter fields, the moduli fields can be dynamically stable. Away from the minimum, 
the evolution of the moduli fields can be expressed in terms of an effective scalar field having 
a double exponential(exponent of an exponential) potential as considered in eq. (3.1) [8]. 

In the next section, we introduce the notations and conventions. Further, keeping 
in sight the ensuing cosmological dynamics, we outline characteristic features of a scalar 
held. Thereafter cosmological evolution equations have been non-dimensionalized for the 
convenience in the dynamical analysis. Critical points for the resulting autonomous system 
and their stability issues are discussed in Sec. Ill, followed by results and conclusions in Sec. 
IV. 

2 Dynamics of canonical scalar field 

The action of a minimally coupled scalar held in a four dimensional spacetime is given by 

S = J d A Xy/^ [ A + ^ + C m ] , (2.1) 
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where g is the determinant of the metric, R is the Ricci scalar, k 2 = 87 tG (c = 1) , C m is 
the matter Lagrangian, and the scalar field Lagrangian £<* is given by 

^ = ( 2 . 2 ) 

with V being a general potential for 4>. Epsilon e is equal to +1 for the quintessence, and — 1 
for the phantom field [9]. The variation of the metric gives the gravitational field equations: 

G\iv = + r ^fj.v[4>)) (2-3) 


where = R^ v — \g^ u R is the Einstein tensor. The matter stress-energy tensor for a 
fluid may be parametrized in terms of the four velocity of the fluid u^, and the density and 
pressure functions p m and p m as 

T//z/ \rn\ {Pm + Pm)u fl Uv+pg flu . (2.4) 

The scalar field energy-momentum tensor is 

= td^d v <\> + . (2.5) 


We consider a spatially flat FRW metric with the metric tensor given by 

9txv = diag(— 1 , a(t) 2 ,a(t) 2 , a{t) 2 ) ( 2 . 6 ) 


where a{t) is the scale factor. The scalar held is taken to be spatially homogeneous, i.e., 
4> = This gives 

< w = -n = -\i 2 -v(<l>), (2.7) 

T m = t a*\ = T m = r* = - W) ■ (2.8) 

We assume that generic equation of state for the ith component (radiation, dark matter or 
scalar held) is 

Pi = (rt~ 1 )Pi (2-9) 


where 7 is a constant. Gamma ( 7 ) is equal to 4/3 for radiation, and 1 for pressureless dark 
matter. For the scalar held, relation (2.9) gives the effective equation of state as 


p<j) + P4> _ 2 <j) 2 

P<t> </> 2 + 2V ' 


( 2 . 10 ) 


Eq. (2.3) relates the Hubble parameter H = d/a to the stress-energy tensor components to 
give the Friedmann constraint 


H 2 


2 (Pm T P(j> T Prad) j 


( 2 . 11 ) 
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and the Raychaudhuri equation 


• K? 4 

H = — — {Pm + 2 Prad + Z<p) ■ 

(2.12) 

Finally, the conservation equations for the field (<)>), cold dark matter(] 
read 

m) and radiation(rad) 

Pm 4 “ 377p m — 0 , 

(2.13) 

p<t> + 377/9^70 = 0 , 

(2.14) 

Prad 4 “ 477 Prad — 0 . 

(2.15) 


So far we have not considered any coupling between the dark matter and the scalar 
field. Introducing an effective coupling denoted by Qj between dark energy and dark matter 
modifies Eqs. (2.13) and (2.14) to 


Pm 4“ p m — Ql > 
P<t> + = -Q1 


while the modified Klein-Gordon equation is given by 


4> + 3 Hcj) + 


dV 

d<j) 


Qi 

e 4> 


(2.16) 

(2.17) 


(2.18) 


Eqs. (2.12), (2.16), (2.17) and (2.18) are the evolution equations subject to the Friedmann 
constraint (2.11). Consider dimensionless variables [5, 8, 10]: 


_ K(j> _ kVV _ Hy/PracL _ — 1 dV 

~ VQH' V ~ VZH' “~ y/3H' ^ ^ VK d(j) ' 


(2.19) 


In terms of these new dimensionless variables, the evolution equations lead to a general¬ 
ized four dimensional autonomous system 


~e + ^j\ y2z + ~ 6x2 ~ y2 ~ ^ + + 2ex2 ^ ~ ' 

( 2 . 20 a) 

R 3 4 

J 7 jV xz + -jVid - £X 2 -y 2 - u 2 ) + -u 2 + 2ex 2 }. 

( 2 . 20 b) 

3 4 

2 u 4 - -u[(l — ex 2 — y 2 — u 2 ) 4 - - u 2 4 - 2 ex 2 ]. 

Z O 

( 2 . 20 c) 

d6z 2 x(r — 1 ). 

( 2 . 20 d) 
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Here prime denotes differentiation with respect to the new time variable N = ln(a/ao) and 
we have further defined 


1 dV 
kV dcj) 


( 2 . 21 ) 


T = V 


d?V 

dip 2 


(dV\2 ’ 
v dtp > 


( 2 . 22 ) 


^ _ k Qi 
** ~ ‘ 


(2.23) 


As the dynamical system (2.20) is invariant under y —» —y, we need to solve the system only 
for y > 0. 

The Friedmann constraint (2.11) gives 


n 2 pm 
3 H 2 


i 2 , 2,2 -1 

+ ex +y +u = 1 


(2.24) 


or 


where 


+ ^777, — 1 




= ft 2 ;00 

~ 3 H 2 ’ 




K Prad 

3 H 2 ’ 


_ K Pm 

= 3H 2 ' 


(2.25) 


(2.26) 


Also is bounded, i.e., 0 < ex 2 + y 2 < 1 for non-negative H r and H m . 

The dynamical system framework (2.20) introduced so far is general enough to accom¬ 
modate analysis for large class of potentials with or without any interaction. Furthermore, 
the nature of the field (whether quintessence or phantom) has been properly incorporated 
into the framework. In the next section we shall work with a specific kind of potential and 
interaction term. Once we know the form of interaction Qj and the potential V((j>), we can 
obtain the critical points (x*, y*, u*, z*) of the autonomous system ( 2 . 20 ) by imposing the 
conditions X* = = u* = = 0. These critical points must be real to ensure their exis¬ 

tence in the phase space. To study the stability of the critical points, let us consider small 
perturbations 5x,5y,8u and 8z around the critical points (x*, y*, it*, z*), i.e., 

x —^ x* T 5x, y —> y* + 8y, u —» u* + 5u, z —>■ z* + Sz. (2.27) 

Combining the four variables into a vector x a , we have x a —)• x a * + dx a . On substituting 
perturbed variables from (2.27) into the system (2.20), we get a set of first order differential 
equations symbolically represented by 

-4(5x a ) = M Xa ,(<5x a ), (2.28) 
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Figure 1. (Color online). Dynamics of scalar field in the phase space when there is no interaction 
for a double exponential potential for different initial conditions. Convergence of both curves is at the 
same point, (c) Evolution of different energy densities f \ w.r.t. N and (d) evolution of w e ff w.r.t. 
N. The initial conditions for the preceding evolution were Xi = 0.015, y; = 0.015, Ui = -\/0.85 and 
Zi = 40, and A = 6. 


where <5x a and M Xa „ are the column vector of the perturbations and the coefficient matrix, 
respectively. M Xat depends on the critical points and is represented as 


dx' 

dx' 

dx' 

dx' 

\ 

dx 

dy 

du 

dz 


d£ 

d</_ 

dyi 

dy' 


dx 

dy 

du 

dz 


du' 

du' 

du' 

du' 


dx 

dy 

du 

dz 


dz' 

dz' 

dz' 

dz' 


dx 

dy 

du 

dz 

/ 


The nature of the four eigenvalues of the coefficient matrixM Xat determines the stability of 
the critical points. The criteria for establishing the stability of the critical points is as follows: 

а) S node : if all the eigenvalues are negative. 

б ) Unstable node : if all the eigenvalues are positive. 

c)Saddle point : if some eigenvalues are positive and some are negative. 


However, if any one or more eigenvalues vanish, the linear stability analysis breaks down. 
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Figure 2. (Color online). Dynamics of scalar field in the phase space when interaction is switched 
on for a double exponential potential for the same initial conditions as in figure 1. Convergence of 
both curves is at the same point. Also as in figure 1, the evolution of different energy densities (c) 
f 1i, and (d) w e ff w.r.t. N, respectively. 


In that case, the center manifold theory needs to be applied [11, 12]. The stability for the 
special case can also be determined from the signs of other non-null eigenvalues if the critical 
points belong to normally hyperbolic sets [13]. 


3 Case of double exponential potential 


Upto now, we have not considered any specific form of V and Qi. Now onwards we shall work 
with e = 1 (quintessence) for the sake of simplicity. Moreover, we shall assume that Qj is 

•y/|which is also a simple and widely studied interaction in a variety of potentials [4, 
5, 9], 


We choose the double exponential potential for our quintessence field. Although this 
potential is significant from the viewpoint of fundamental theories, particularly of higher 
dimensions [7], it is still under-studied potential in the literature for dark-energy modeling. 
It is to be mentioned that this potential is different from the sum of two exponential terms [13] 
which is also known by the same name. 
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On putting potential 


V(<f>) = V 0 exp(e~ KX< f > ) 


(3.1) 


in Eqs. (2.21)-(2.22), we get 


z = A exp 


—k\( j) 

A 


r -1 = - 

z 


Also, using Qi = \/ 3 


in eq. (2.23), we get 

Q = f3 (1-x 2 -y 2 - u 2 ) 


(3.2) 

(3.3) 

(3.4) 


Therefore, the dynamical system (2.20) transforms as 

x' = -3x + \ j^y 2 z + |x[(l - x 2 - y 2 - it 2 ) + ~u 2 + 2x 2 ] - f5 (1 - x 2 - y 2 - u 2 ) (3.5) 

V = -\J^yxz + |y[( 1 - X 2 - y 2 - u 2 ) + ^u 2 + 2x 2 ] (3.6) 

u' = —2 u + -u[(l — x 2 — y 2 — u 2 ) + -u 2 + 2x 2 ] (3.7) 

2 3 

z = ~V6zx\ (3-8) 

Also, from Eqs. (2.11) and ( 2 . 12 ), we have 

^ = -|[l +»«//] (3-9) 

where w e ff in terms of dimensionless variables (2.19) is 

w e ff = 2x 2 — 1 + -u 2 + (1 — x 2 — y 2 — u 2 ) . (3.10) 

3 

To have acceleration, we need w e ff < — 

The auxiliary variables a, b, c, /rl, /r2, /r3 and ^4 in table 2 are defined as 


/xl = ^C9 2 + >/3%/^9 i ), 

A*2 = ^(-/3 2 + V3V/3 2 - /? 4 ), 

^3 = ^2 [/ 5 2 + VP 2 (9 + 5 / 5 2 )], 

/x4 = /3 2 + /3 2 (9 + 5/3 2 )], 

4/1 2 + 9 

a= g > 

, 4/3 2 - 9 

4/3 2 - 3 


(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 
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Table 1. Critical points of the system and their relevant properties for double exponential potential. 


pt. 

x c 

Uc 

Uc 

Zc 


7* 

U<f, 

Existence 

A 

0 

0 

1 

Z 

0 

Indet. 

Indet. 

V/3 

B 

-1 

0 

0 

0 

1 

2 

1 

V/3 

C 

0 

0 

1 

0 

0 

Indet. 

Indet. 

V/3 

D 

1 

0 

0 

0 

1 

2 

1 

V/3 

E 

i 

2 0 

0 

V-3+4/3 2 

2/3 

0 

1 

w 

2 

1 

/3 2 >f 

F 

2/3 

3 

0 

0 

0 

4/3 2 

9 

2 

1 

V/3 

G 

0 

1 

0 

0 

1 

0 

-1 

V/3 

H 

3 

v/i/P+9 

0 

0 

1 + A 

2/3^ 

9 

—2/3 2 

/3 2 >^ 

2/3 

2/3 

2/3 2 +9 

2/32+9 


Table 2. Eigenvalues for the critical points, their stability properties and energy densities of various 
components. See text for a, b, c, pi and p2. 


pt. 

“»«// 

q 

S 

i-fG 

E.V’s 

Hr 

r 2 m 

Stability 

A 

l 

3 

l 

Indet. 

( 2 ,- 1 , 0 , 1 ) 

1 

0 

Saddle 

B 

1 

2 

Indet. 

(3,1,3 + 2/3, V6A) 

0 

0 

Unstable 

C 

1 

3 

l 

Indet. 

( 2 ,- 1 , 0 , 1 ) 

1 

0 

Saddle 

D 

1 

2 

Indet. 

(3,1,3 — 2/3, —\/ 6 A) 

0 

0 

Saddle 

E 

1 

3 

l 

1 

2 

(pl,p2,2,~^) 

1 - A 

1 4^ 

1 

w 

Saddle 

F 

4/3 2 

3+4/3 2 

4/3 2 

(a, b, c, —2y^|/3A) 

0 

1 4/3 2 

9 

Saddle 


9 

6 

—4/3 2 +9 

G 

-1 

-1 

Indet. 

(-3,-2, 0,-3) 

0 

0 

stable 

H 

-1 

-1 

1 2/3 2 

9 

(A/i4,-2,^) 

0 

-9 

w 

Saddle 


Now we shall scrutinize each critical point carefully: 

a) From table 1, it is clear that the constraint /3 2 > ^ for the existence of point H can 
not be satisfied with the constraint equation which needs < 0 for any real /3. Therefore 
this point is excluded from further consideration. 
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b ) As w e ff = | and <7 = 1 for points A, C and E (see table 2), the scalar field mimics the 
nature of radiation while the universe is decelerating for these points. As these are the saddle 
points too, this is just a transient behavior. 

c) As = 1 and q = 1 for unstable points B and D (see table 2), these are the dark energy 
dominated solutions with the universe decelerating. 

d) The only attractor we have in our case is the point G. This point has w e ff = — 1 (from 
figure 1(d)) and q = — 1, which are the conditions to be satisfied for accelerated expansion 
(see eq. (3.10)). On putting w e ff = —1 in eq. (3.9), H = 0 vanishes, in other words, H 
becomes constant, which leads to a oc e Ht . This implies that we have a de-Sitter solution at 
the present epoch. No scaling solution is observed for this point since f= 1. However, the 
interesting aspect of this point is that we have a global attractor at the present epoch. It 
means that our trajectories converge at this point independent of the initial conditions (see 
figure 1(a) and figure 1(b)). 

From figure 1(c), it is apparent that when there is no interaction (i.e., /3 = 0), the above 
defined autonomous system follows the desired cosmological dynamics. It describes an early 
radiation-dominated era(RDE, dominated by £l ra d) followed by a matter dominated epoch 
(MDE,dominated by Il m ), and finally the universe enters into the present accelerated regime 
(dominated by Q^). 

4 Conclusion and outlook 

In this work, the deSitter solution is obtained at the present epoch to account for the ac¬ 
celerated regime. Dynamical evolution of the field translates into a smooth transition of the 
universe from a decelerated epoch to an accelerated one in the presence of a dark energy 
component — a scalar field having a self interaction described in terms of a double exponen¬ 
tial potential. Suitable combinations of initial conditions and parameter values are found for 
which a realistic sequential dominance of different cosmological constituents could follow. It 
is found that the initial dynamical evolution is dependent on the interaction parameter (5; 
however, the asymptotic behavior of the system is independent of the value of (3 — typical 
of an attractor solution(see figure 1(c) and 2(c)). It is so because there exists only one stable 
attractor, i.e., point G. As the point G admits a deSitter solution (complete DE dominated), 
it is not possible to address the coincidence problem strictly within the scenario considered. 

As alternative direction for future work, it would be interesting to explore if other 
possibilities such as time-dependent interaction, non-minimal coupling and phantom field 
(usually considered in the literature [4, 5, 9]), might provide a global attractor for which HjS 
are parameter dependent and which could account for the coincidence problem. 

5 Acknowledgement 

We are grateful to several members of Department of Physics and Astrophysics for their useful 
comments. We would also like to thank Shiv Sethi for a critical review of the manuscript, 
and Sanjeev Kumar for help in numerical methods. VG wishes to thank the CSIR (India) 
for the financial support through grant number 09/045/(0933)/2010 EMR-I. 

References 

[1] S. Perlmutter, G. Aldering, G. Goldhaber, R. A. Knop, P. Nugent, P. G. Castro, S. Deustua, 

S. Fabbro, A. Goobar, D. E. Groom, I. M. Hook, A. G. Kim, M. Y. Kim, J. C. Lee, N. J. Nunes, 


- 9 - 


R. Pain, C. R. Pennypacker, R. Quimby, C. Lidman, R. S. Ellis, M. Irwin, R. G. McMahon, 

P. Ruiz-Lapuente, N. Walton, B. Schaefer, B. J. Boyle, A. V. Filippenko, T. Matheson, A. S. 
Fruchter, N. Panagia, H. J. M. Newberg, W. J. Couch, and T. S. C. Project, Measurements of 
?? and ?? from 42 high-redshift supernovae, The Astrophysical Journal 517 (June, 1999) 565. 

[2] V. Sahni and A. Starobinsky, The case for a positive cosmological lambda-term, International 
Journal of Modern Physics D 09 (Aug., 2000) 373-443. 

[3] V. Sahni, Dark matter and dark energy, arXiv:astro-ph/0403324 (Mar., 2004). 

[4] E. J. Copeland, M. Sami, and S. Tsujikawa, Dynamics of dark energy, International Journal of 
Modern Physics D 15 (Nov., 2006) 1753-1935. 

[5] C. G. Boehmer, G. Caldera-Cabral, R. Lazkoz, and R. Maartens, Dynamics of dark energy with 
a coupling to dark matter, Physical Review D 78 (July, 2008). 

[6] E. J. Copeland, A. R. Liddle, and D. Wands, Exponential potentials and cosmological scaling 
solutions, Physical Review D 57 (Apr., 1998) 4686-4690. 

[7] G. Huey, P. J. Steinhardt, B. A. Ovrut, and D. Waldram, A cosmological mechanism for 
stabilizing moduli, Physics Letters B 476 (Mar., 2000) 379-386. 

[8] S. C. C. Ng, N. J. Nunes, and F. Rosati, Applications of scalar attractor solutions to cosmology, 
Physical Review D 64 (Sept., 2001). 

[9] X.-m. Chen, Y. Gong, and E. N. Saridakis, Phase-space analysis of interacting phantom 
cosmology, Journal of Cosmology and Astroparticle Physics 2009 (Apr., 2009) 001-001. 

[10] L. Amendola, Coupled quintessence, Physical Review D 62 (July, 2000). 

[11] A. A. Coley, Dynamical Systems and Cosmology. Springer Science & Business Media, Oct., 
2003. 

[12] J. Wainwright and G. F. R. Ellis, Dynamical Systems in Cosmology. Cambridge University 
Press, June, 2005. 

[13] G. Leon, Y. Leyva, and J. Socorro, Quintom phase-space: Beyond the exponential potential, 
Physics Letters B 732 (May, 2014) 285-297. 


- 10 - 



